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Abstract 

We analyze a model of hypercubic random surfaces with an extrin- 
^P ' sic curvature term in the action. We find a first order phase transition 

at finite coupling separating a branched polymer from a stable flat 
phase. 



1,. 1 Introduction 

(D ■ 
,^ ■ Random surfaces have attracted a lot of interest in different branches of 

physics in the recent years. The thermodynamics of surfaces embedded in 
"real" three-dimensional Euclidean space is interesting in itself as it may 
describe some properties of real membranes appearing in nature. In high 
5r , energy physics random surfaces play a role for example as the world-sheet 

of strings or as the space-time of two-dimensional quantum gravity. In this 
paper we investigate the hyper-cubic random surface model with an extrinsic 
curvature term in the action. Although originally motivated by bosonic 
string theory it may well be interpreted as a fluid phantom membrane with 
bending rigidity. The word phantom refers to the fact that the surface 
considered here can, differently form real- world surfaces, self- intersect. 

The hypercubic random surface model was originally proposed by Wein- 
garten Q as a non-perturbative regularization of the world-sheet of bosonic 
strings. It was, however, soon realized that the dominating geometries in 
this model have the essentially one dimensional structure of branched poly- 
mers. It was shown M that for spherical surfaces the string susceptibility 



exponent 7 is equal to 1/2, the generic value for branched polymers. In |^] 
the model was extended to include an extrinsic curvature term and it was 
shown, under some assumptions, such a term in the action does not change 
the model's critical behavior for any finite coupling. It therefore came as a 
surprise when numerical evidence for a non-trivial behavior with 7 = 1/4 
was observed [Q when an additional local constraint on the allowed configu- 
rations was introduced, a term forbidding self-bending surfaces. The result 
7 = 1/4 was especially interesting because this value fits nicely into a series 
7 = 1/n, n = 2, • • • of positive string susceptibility exponents discussed in 
[^. Hence this result led to some speculations Q. However, using improved 
numerical techniques and larger lattice-sizes it was shown ^ that the true 
large volume behavior is masked by strong finite size effects but nonetheless 
is well described by branched polymers with 7 = 1/2. 

In this paper we generalize the hyper-cubic random surface model by 
adding a term coupled to the external curvature to the action. Introducing 
such a term will certainly change the dominating geometry, at least in the 
infinite coupling limit. In this limit the dominating geometry is flat with 
external Haussdorff dimension^ dh = 2 different from dh = A for branched 
polymers. We will demonstrate numerically that the transition actually 
happens at finite coupling. This means we observe for a fluid membrane 
a crumpling transition separating the branched polymer from a stable flat 
phase. 

On the first glance, the announced transition at finite coupling seems to 
be ruled out by 0, which states that "the critical exponents take their mean 
field value if the susceptibility of the model and a coarse grained version of 
it both diverge". In other words under the assumptions 

7 > and 7 > (1) 

the model is always in the branched polymer phase for finite couplings. 
This result was obtained with a renormalization group argument where the 
branched polymer was decomposed into " blobs" , components which can not 
be cut into two parts along a loop of length two. A key step in the derivation 
of this result in [^ is the formula 

;,M = .^(a,^, (2) 

which relates the susceptibility x{^^) of the original model to the susceptibil- 
ity x(A) of the decomposed blobs with renormalized coupling Jl. Obviously 
^for a definition see eq. [nl 



X = 1 if the susceptibility x diverges. Under the assumptions stated above, 
i.e. 7 > for the blobs one concludes that the blobs are not critical p. > fie 
and X is analytic at this point. With a Taylor expansion one gets the self 
consistency relation 7 = 1 — 7 which is solved by 7 = 2 > the generic value for 
branched polymers. However, Durhuus emphasized [^ that the condition 
X = 1 alone does not imply that the blobs are non-critical. If 7 < the 
susceptibility x does not diverge and x = 1 can be satisfied at the critical 
point fi = fie oi the coarse grained system. At such a point the entropy is 
still dominated by the branching of the surface but the branches themselves 
are critical, which in effect changes the exponent 7 for the whole system. If 
one assumes that 7 takes the KPZ-values [^] one can derive from (g) a series 
of positive 7 = 1/n, n = 2, • • • |5|. 

The derivation of this series is rather formal. It does not provide a 
description how to construct a system which has this property. We want 
to check numerically if the hyper-cubic random surface model exhibits non 
trivial behavior at an eventual crumpling transition. In other words we check 
if the assumptions (|l|) are fulfilled in the whole range of couplings eg. Non- 
trivial behavior can be expected only if these assumptions are not satisfied 
at a possible phase transition induced by the external curvature coupling. 
A special focus will therefore be given to a possible transition point. 

2 The model 

A hyper-cubic surface is an orientable surface embedded in Z^ , obtained by 
gluing pairwise together plaquettes P along their links until no free edges are 
left. A plaquette P is a unit-square which occupies one of the unit-squares 
'S'zM, /U = 1- ■ ■ D in the embedding lattice. In the following we use the word 
square to refer to the embedding lattice and the words link or plaquette to 
describe the internal connectivity of the surface. A link is shared by exactly 
two plaquettes, which are glued together along each one of their edges. Note 
that a link or a square in the embedding lattice can be occupied more than 
once which means these are non self-avoiding, i.e. phantom, surfaces. 
The partition function is 

^(/3) = ^e-^ = ^e-^Ar(A), (3) 

EeS A 

where the sum runs over an ensemble S of hyper-cubic surfaces with fixed 
(spherical) topology. The number A/'(A) of surfaces with a given area A is 

^see eq. (pf) 



expected to grow like 

M{A) ^ ef'^^A^-^ (4) 

in the large volume limit. The coefficient 7 is the entropy exponent. To 
balance the exponential growth of the number of surfaces with the volume, 
a term ^A linear in A has to appear in the action. The partition function 
is well-defined only if ^ < ^c- For the extrinsic curvature we note that 
there are three possible ways to embed two neighboring plaquettes pi,P2 in 
a hyper-cubic lattice Z^. The external angle G^, which we assign to link /, 
can therefore take three possible values: 

0: pi and p2 occupy two neighboring squares in the same coordinate 
plane, 

^: pi and p2 occupy two neighbor squares in diff'erent coordinate planes, 

tt: Pi and p2 occupy the same square. We call such a configuration self- 
bending. 

We use the potential 

Se.t[E]ie) = 11A+ ^Y,{1 - cosQi) = 11A + eY, [^n,e, + 1^/2,0^ ■ (5) 

The physical questions we want to address are most conveniently ana- 
lyzed using the canonical ensemble. However, the algorithm used for this 
work Q requires moderate volume fluctuations for ergodicity. Therefore we 
simulate the quasi-canonical ensemble: 

Z(A,e) = Y, e-^[^l(^'^) = Y <5^,^e-^[^l(^'^)-^(^-^)' (6) 

Ee&^ Ees 

In the first expression, defining the canonical ensemble, the sum runs over all 
surfaces with area A. In the actual simulation surfaces with areas different 
from A are generated, where the area-fluctuations are conflned by the addi- 
tional Gaussian potential. The information about the canonical ensemble is 
extracted by measuring only if A = A. 

2.1 Numerical results 

To localize possible transitions we measure the energy fluctuation of the 
external curvature field 

^ d\ ^, , , < G2 > - < e >2 

Xe.t = ]4 ^ log Z{A, e) = J (7) 
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Figure 1: The susceptibility Xext in D = 3 dimensions 
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Figure 2: The distribution of the average external curvature @ in D = 3 
dimensions obtained with re-weighting. The non zero probability for un- 
physical negative values G) < is an artifact of the extrapolation to this 
region in the re- weighting process. 
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Figure 3: The pseudo-critical coupling ecrit{^/A). 



and search for maxima in this observable. The average external curvature 
is proportional to the first derivative of the free energy : 



<e>=i^|'°s2(A,.) 



^EEe-»-'^'. 



2A 



(8) 



We simulated in three and four embedding dimensions and scanned the 
coupling range < e < 3 for maxima of Xext- In both cases we found a 



single peak. We used re-weighting methods |l^ to extract the shape of the 
peaks shown in figure ^ from four independent measurements per volume at 



e ^ epc in d = 3. 



The peaks grow quickly, in fact a bit faster than linear, with the volume. 
This suggests a first order phase transition, which is confirmed by a look at 
the distribution of external curvature 0, where the bar indicates the average 
taken over a given lattice. In figure we show this distribution for D = 3 



and the coupling e 



'^pc 



close to the pseudo-critical coupling. One observes 



a clear signal of a first order phase transition, namely two maxima separated 
by a minimum which becomes deeper as the size of the surfaces is increased. 
This clearly indicates two separate phases. In one phase the average external 
curvature is close to zero, the typical surface is flat. 

To estimate the critical coupling we use that for a first order phase transition 
one can expect 

ej,,{A)=e, + j + 0{j^) (9) 

for the scaling of the pseudo-critical couplings Cpc- In figure y we show the 
numerical estimates for e„c- Note that the pseudo critical coupling decreases 
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Figure 4: The schematic phase diagram for the Hyper-cubic random surface 
model with extrinsic curvature. One finds qualitatively the same behavior 
in D = 3 and D = A embedding dimensions. 

with the volume which indicates that Ec is finite in the thermodynamic limit. 
With a linear fit to equation (|9|) we find for the critical coupling 

D = 3 ec = 1.530(3) 
D = A ec = 1.643(4) ' 

To demonstrate the change in the geometrical behavior we consider the 
radius of gyration 

G\A, e) =< i ^{S,{p) - S,ipo)f >, (10) 



A 



p 



which is the average squared distance of the plaquettes to some reference 
plaquette po. The S^ are the coordinates of the plquette in the embedding 
space. In the numerical simulations we do a stochastic average over the pq 
by choosing 10 % of the plaquettes at random as po pei' measurement. One 
can use the radius of gyration to define the external Haussdorff dimension 

dh 

G'^{A,e)^A^/'^'^. (11) 

In a sense the Haussdorff dimension is the largest possible dimension for the 
embedding space which can still be completely filled by the surface. We have 



measured numerically the radius of gyration for e = and e = 2 for volumes 
A in the range 200 to 5000. With a fit to (pA|) we extract the Haussdorff 

dimension: 

e = 0. dh = 4.07(3) 
e = 2. dh = 2.00(3). 

The value found for e = is two standard deviations away from the expected 
result Q but this should presumably be attributed to finite size effects. 

The results are summarized in the phase diagram figure Q. The model 
is defined only for fi > /ic(e) which defines the critical line depicted by the 
solid line in the phase diagram. Above this line we find two phases. For 
e < ^crit the system is in the branched polymer phase with df^ ~ 4. If we 
assume hyperscaling 

u = dj,\ 7 = 1.(2-,?) (12) 

the anomalous scaling dimension is ry = as expected for a branched poly- 
mer. The strong coupling phase is flat, which is reflected by the observation 
dfi = 2. In this phase basically no baby universes are present. For example, 
on a surface of size A = 500 we found at e = 2 on average only 0.06 baby 
universes of size B = 5. The largest one found in this phase is of size 33 
and appeared only once in 200000 measurements. For comparison in the 
branched polymer phase at e = we found in average 10.6 baby universes 
of size 5 on the surface, the probability for the largest possible one with 
B = 249 is about 0.02. Therefore one cannot use the baby universe distri- 
bution to measure 7 in the flat phase. However the fact that there effectively 
are no baby surfaces means that 7 is large negative. Therefore the existence 
of this phase is not in contradiction with |Q because the assumptions (|l|) 
are not met. 

3 Summary 

We have analyzed the model of hypercubic random surfaces with an extrinsic 
curvature term in the action. We observed a first order phase transition, 
which separates a branched polymer phase from a flat phase. We have shown 
that the critical coupling is flnite in the thermodynamic limit, i.e. the flat 
phase survives in the large volume limit and we observe true long-range 
order. This may seem a bit surprising flrst because long range correlations 
in two dimensional systems are rather unusual. In a comparable model of 
triangulated fluid membranes the existence of a stable, long range ordered 



phase had somewhat been disputed |11]. 



Our hope to find a non-trivial, positive 7 for the geometry at the transi- 
tion of the external-curvature "field" was not satisfied. Although we find a 
phase transition, we do not have the situation where the exponent 7 of the 
individual branches is negative while the exponent for the overall geometry 
is positive but different from the generic branched polymer value 7=5- 
Instead we find 7 large negative, which is also the reason why the transition 
is not in contradiction with the statement "no transition at finite coupling" 
for this model given in H, because this argument relies on the assumption, 
that 7 is positive. 
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